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Physical phenomena driven by topological properties, such as the quantum Hall effect, have the 
appealing feature to be robust with respect to external perturbations. Lately, a new class of materials 
has emerged manifesting their topological properties at room temperature and without the need of 
external magnetic fields. These topological insulators are band insulators with large spin-orbit 
interactions and exhibit the quantum spin-Hall (QSH) effect. Here we investigate the transition 
between QSH and normal insulating phases under topological deformations of a two-dimensional 
lattice. We demonstrate that the QSH phase present in the honeycomb lattice looses its robustness 
as the occupancy of extra lattice sites is allowed. Furthermore, we propose a method for verifying our 
predictions with fermionic cold atoms in optical lattices. In this context, the spin-orbit interaction 
is engineered via an appropriate synthetic gauge field. 
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I. INTRODUCTION 

In the last three decades the physics community has 
been fascinated by topological states of quantum mat- 
ter, initiated by the discovery of the quantum Hall ef- 
fect (QHE) by von Klitzing in 1981 [1]. Recently it has 
been discovered that a new class of materials, referred 
to as topological insulators, show robust conducting edge 
states even in the absence of external magnetic fields [2- 
4]. These materials are two- or three-dimensional band 
insulators with a large spin-orbit interaction (SOI) and 
exhibit the so-called quantum spin-Hall (QSH) effect [5- 
8]. Here, spin-orbit effects play a role similar to the ex- 
ternal magnetic field in the QHE. However, contrary to 
an external magnetic field, SOIs conserve time-reversal 
symmetry (TRS), implying that the edge states appear as 
Kramers doublets or equivalently as counter-propagating 
modes along the system boundary carrying opposite spin. 
In two-dimensional (2D) systems a Z2 topological invari- 
ant v specifies the robustness of these helical edge-states 
in the presence of disorder and interactions [9]. While 
v = 1 states a topologically protected (i.e. stabilized) 
QSH phase, v = indicates that the edge states are un- 
stable with respect to impurity scattering therefore re- 
ducing the system to a normal insulator (NI). 

In this Article we explore the interplay between spa- 
tial lattice topology, which is characterized by the lattice 
connectivity, and the topological order associated to the 
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FIG. 1: (Color online), (a) HCL with lattice vectors vi and 
v 2 and two atoms in the unit cell, (b) 75 lattice with three 
atoms in the unit cell. Rims A and B (open circle and square) 
have a lower connectivity of 3 compared to 6 for the hub site 
(filled circle). The 75 lattice can be viewed as two nested 
HCLs (emphasized by the red and blue hexagon). 



quantum spin-Hall phases. The latter is characterized 
by the Z 2 topological invariant defined on a fibre bundle 
associated to the wave- functions over the first Brillouin 
zone [2]. We show that spatial deformations — modi- 
fying the lattice topology — induce a change in the Z2 
invariant, thus leading to a topological phase transition 
from a QSH insulator to a NI. 

Specifically, we study the effect of adding an additional 
site to the unit cell of a honeycomb lattice (HCL) sub- 
jected to SOI. While this lattice deformation preserves 
the relativist ic behavior of the low-energy modes, we 
show that the system undergoes a topological phase tran- 
sition that changes its topological class, labeled by the Z2 



index, and its transport properties. The possibility to in- 
vestigate this topological phase transition with ultra-cold 
atoms trapped in optical lattices is promising [10, 11]. In 
this context, cold- atom experiments are not only suitable 
for engineering the QSH phase [12, 13] and to test its sta- 
bility with respect to controllable disorder [14], but also 
to induce lattice transformations. 



II. MODEL AND FORMALISM 



Adding an extra site to the unit cell of a HCL gives 
rise to the so called % lattice, c.f. Fig. 1. The 75 lattice 
not only differs from the HCL by the number of inequiva- 
lent lattice sites but also by their respective coordination 
numbers: while the HCL (Fig. la) is characterized by 
two inequivalent lattice sites A and B, each with coordi- 
nation number 3, the Ts lattice (Fig. lb) has a unit cell 
with three inequivalent sites. Two of the latter (generally 
called rim A and B) are characterized by a coordination 
number 3 while one site (referred to as hub H) is con- 
nected to 6 nearest-neighbors [15-17]. 

Both lattices are modeled by a tight-binding (TB) 
Hamiltonian Ho = t^2( n ^m)a c lia c ma with spin inde- 
pendent nearest-neighbor hopping amplitude t. Here, 
c lia( c na) is the creation (annihilation) operator for a par- 
ticle with spin direction a on the lattice site n. Regarding 
the SOI, symmetries allow for two types of interactions 
that preserve TRS, namely the Rashba-type [2] and the 
Haldane-type [18]. The effects of the Rashba-type SOI on 
the transport properties of the HCL and 75 lattice have 
been reported in Refs. [2, 20]. The Haldane-SOI, which is 
the interaction of interest in the present study, is modeled 
via a spin-dependent second-neighbor hopping term 



Hso = i^so ^2 ^2 c ^« ( di x d i) ' a <*P Cm ^ ' ( 1 ) 

a(3((n,m)) 



Here, the cr^ are matrix elements of the Pauli-matrices 
<y with respect to the final and initial spin states a 
and j3 and d^- are the two displacement vectors of the 
second-neighbor hopping process connecting sites n and 
m. Since in 2D lattices hopping is naturally restricted to 
in-plane processes, the SOI is effectively proportional to 
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FIG. 2: (Color online), (a) Energy spectrum as a function 
of the in-plane momentum (k x ,k y ) and for a fixed value of 
the spin-orbit interaction £so — 0.25£. The first Brillouin 
zone and its six corners are indicated, (b) Cut of the energy 
spectrum at k x = (shadowed plane in (a)). 



III. TRIVIAL VS. TOPOLOGICAL 
INSULATORS 



We start by analyzing the effect of finite SOI on 
the energy spectrum of the infinite T3 lattice. Because 
of the unequal connectivity of rim and hub sites, the 
Hamiltonian (1) effectively induces hopping between rim 
sites of the same kind while this is cancelled for hub 
sites. The spectrum is obtained by exact diagonaliza- 
tion (cf. App. A) and the result is shown in Fig. 2. Due 
to its topological origin the non-dispersive band at the 
charge neutrality point is not affected by SOI. On the 
contrary, two bulk energy gaps open between the non- 
dispersive and the electron/hole branches, respectively. 
This property can be directly deduced from the single- 
valley long- wavelength approximation of the Hamiltonian 
(cf. App. B). We find that the general form of this rela- 
tivists Hamiltonian — describing both the HCL and the 
73 lattice — reads 



U = v F 5] • p <g> I 2 + Aso£ 2 ® cr z 



(2) 



In the absence of SOI the energy spectra of the two 
lattices are characterized by two identical, electron-hole 
symmetric branches [2, 19]. Moreover, for the T3 lattice a 
unique non-dispersive band is present at the charge neu- 
trality point [19]. This is rooted in the lattice topology, 
which allows for insulating states with finite wave function 
amplitudes on the rim sites and vanishing amplitudes on 
the hubs. 



Here v^ is the Fermi velocity, p = — i(d x ,d y ,0) the 
momentum operator (h = 1) in the lattice x?/-plane, 
Aso is the spin-orbit coupling constant, and I2 is the 
two-dimensional identity matrix. The pseudo-spin op- 
erator XI for the HCL is given by the Pauli matrices 
X) = (<Xe, <j y , <j z ), thereby describing a S = 1/2 pseudo- 
spin system. On the contrary, for 75 the pseudo-spin 




FIG. 3: (Color online). One-dimensional energy bands for 
a strip of (a) HCL and (b) 75 lattice with tso = 0.05 1. 
The bands crossing the gap (blue lines) are spin filtered edge 
states. 
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TABLE I: (Color online). In the table the valley-spin is 
associated with the topology of the Bravais lattice and the 
pseudo-spin with the number of lattice sites in the unit cell. 



operator T, = (E x , Y> y , T> z ) is given by [19] 
/n i 




(3) 



These matrices satisfy angular momentum commutation 
relations and describe an S = 1 pseudo-spin (cf. App. B). 

A direct way of visualizing the SOI effects is to consider 
a finite piece of the lattice. The associated TB Hamil- 
tonian can then be diagonalized with periodic boundary 
conditions imposed along one of the spatial directions. 
The HCL with finite SOI is characterized by the QSH 
phase: for each energy value within the bulk energy- 
gap exists a single time-reversed pair of eigenstates on 
each edge of the lattice (see Fig. 3a) — these form two 
Kramers doublets, one on each edge. The conservation 
of TRS prevents the mixing of this couple of states due 
to small external perturbations and scattering from dis- 
order [2]. The QSH phase is intrinsic to the bulk energy 
spectrum and can be validated by topological constants. 
Two important topological invariants are the Z2 index 
v [2] and the spin Chern number n a [21, 22]. In the 
absence of spin- mixing perturbations, both are related 
by v — n cr mod2, where n a = (N^ — N^)/2 and N^^ 
represent the Chern numbers associated to the individ- 
ual spins. For the HCL it has been shown that v — 1 
therefore classifying it as a topological insulator. On the 
other hand it was found that SOI of Rashba-type does 
not induce this kind of edge-states [2] . 

The physical picture is more involved in the case of 
a T3 lattice with SOI. Here gapless edge states are also 
found for energies located within the gaps between the 
electron/hole band and the non-dispersive band, respec- 
tively (see Fig. 3b). However, contrary to the HCL case, 
we find two distinct crossing points for these edge states. 
For each energy value inside the bulk gaps two couples of 
time-reversed Kramers doublets are present, which there- 
fore give rise to four pairs of states along the two edges. 
As a result, even if SOI allow for edge states within the 



bulk energy gap, these are not protected against disor- 
der by TRS — a scattering potential can couple a left- 
mover into a right-mover and vice versa. This is con- 
firmed by the topological invariants (cf. App. C): we find 
that the Z2-index vanishes and that the spin Chern num- 
ber na = 2. For this reason the T3 lattice is topologically 
trivial. 



A transition from the HCL to the T3 lattice is accom- 
plished by varying the hopping parameter between hub 
and rim B from to t. The remarkable result is that 
even for an infinitesimal coupling, the system is in the 
trivial insulating phase. Therefore, we conclude that the 
different behavior of the HCL and the % lattice with re- 
spect to the SOI arises from the different topologies of 
the two systems. The T3 lattice can be viewed as two 
nested HCLs (see Fig. lb). This supports the conjecture 
that two coupled topological insulators behave as a topo- 
logically trivial insulator [22, 23]. We can interpret the 
effect of lattice topology also from a different perspective. 
Although both lattices allow for a description in terms 
of a relativistic Hamiltonian (2) in the long- wavelength 
approximation, only the HCL that is associated with a 
pseudo-spin 1/2 is a topological insulator whereas the 
T3 lattice with pseudo-spin 1 is not. Nevertheless, the 
existence of other topological insulators associated with 
pseudo-spin 1 [24] suggests that the key factor is not the 
dimension of the pseudo-spin but rather the total spin of 
the systems under consideration, cf. Tab. I. In terms of 
the relativistic approximation (2) the total spin is given 
by the combination of the valley-spin, accounting for the 
topology of the Bravais lattice, the pseudo-spin, associ- 
ated with the number of lattice sites in the unit cell, and 
the real fermion spin. In the case of the HCL we have to 
combine three half-integer spins resulting in a half-integer 
total spin. On the contrary, for the T3 lattice we obtain 
an integer spin by combining two half-integer spins with 
an integer one. This conjecture is further supported by 
the fact that both the kagome lattice [25] and the Lieb 
lattice [24] are non-trivial topological insulators, while 
the bilayer graphene is trivial [22, 26]. 
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FIG. 4: (Color online), (a) Phase diagrams in the Fermi 
energy-phase (E-p — 0)-plane for the honeycomb lattice. The 
dark regions represent the energy spectrum. The gaps are 
colored according to the Z2 index: orange [resp. green] gaps 
correspond to v = 1 [resp. v — 0]. (b) Phase diagrams in 
the (Ef — 0)-plane for the 75 lattice, (c) Cut of the phase 
diagram for the HCL in (a) for <j) — 1/21. The orange bars 
denote the energy intervals for which v — 1. (d) Cut of the 
phase diagram for the 75 in (b) for <j) = 1/11. Within the 
relativistic regime limited by i^vHS = 1 [resp. £?vhs ~ 1.4], 
the gauge field (4) opens non-trivial [resp. trivial] Z2 bulk 
gaps for the HCL [resp. 75 lattice]. 



IV. OPTICAL LATTICE VERIFICATION 

Several proposals for realizing optical lattices resem- 
bling the HCL and the 75 lattice have been put forward 
recently. The HCL lattice can be obtained by three copla- 
nar plane waves with the same frequencies and a relative 
angle of 120° [27]. On the other hand, the 75 lattice can 
be realized considering three pairs of counter-propagating 
plane waves with a relative angle of 120° and a relative 
rotation of the polarization plane [28]. Both set-ups share 
the same spatial configuration of the laser beams, there- 
fore an adequate modulation of the laser intensity can 
transform one lattice into the other. The TB Hamilto- 
nian and their long wave length approximations (2) are 
valid for optical lattices populated by fermionic atoms, 
e.g. 40 Kor 6 Li [10]. 

Laser-assisted tunneling is a standard technique for im- 
plementing gauge fields for cold atoms, however a direct 
implementation of the Haldane-type SOI (1) is extremely 
complex. On the other hand, the QSH effect in 2D can 
be interpreted as two standard QHEs acting on each spin 
separately. Therefore, an alternative approach is to de- 
sign an effective vector potential A = (A x ,A y ,A z ) x <j z 
that preserves TRS and opens a bulk energy gap. Re- 
cently, this type of synthetic SU(2) field — practically 
realizable with state of the art techniques — has been 
proposed for a fermionic atom-chip forming a square lat- 



tice [12]. Here we propose to subject the atoms to a 
synthetic SU(2) gauge field 



0,^,0 



X G, 



(4) 



that is associated with an opposite magnetic field for 
spin- up and spin-down electrons. Here S is the area of 
the lattice plaquette and <\> is the number of magnetic 
flux quanta per plaquette felt by each spin, i.e., ±</> for 
spin- up and spin-down, respectively (cf. App. D). Note 
that the gauge field (4) is space-dependent while preserv- 
ing TR symmetry [29]. The TB Hamiltonian including 
this gauge field for both the 75 lattice and HCL in the 
long-wavelength approximation reads 



u = v f j: • p <g> i 2 + a (#) £3, <g> cr z 



(5) 



Contrary to the Hamiltonian (2), the SOI term con- 
tains Yi y and the coupling A(x) = 2k<J)VyS~ 1 x is space- 
dependent. While the latter gives rise to a complex en- 
ergy spectrum as a function of 0, we recover the features 
induced by the Haldane-type SOI in the "low-flux" limit 
(j) <C 1. The Z2-phase diagrams for the HCL and 75 lat- 
tice subject to the synthetic gauge field (4) are shown 
in Figs. 4 (a) and (b), respectively. For arbitrary values 
of </>, the gauge field leads to complex fractal structures 
featuring alternations of non-trivial (orange) and trivial 
(green) phases. However, we note that the HCL has a 
greater tendency to open QSH gaps than the 75 lattice. 
The relativistic regime — described by Eq. (5) — is re- 
covered in the low- flux limit (j) <C 1 [30]. This regime is 
delimited by the van Hove singularities (VHS) ±^vhs i n 
the density of states. In Figs. 4 (c) and (d), we show 
the Z2 index v as a function of the Fermi energy Ep for 
(j) <C 1. In the relativistic regime (yellow area), the HCL 
is characterized by QSH gaps with v — 1, while the 75 lat- 
tice is characterized by trivial insulating gaps with v = 0. 
In this energy range we observe that the synthetic gauge 
field (4) reproduces the effects induced by the Haldane- 
SOI described above. Besides, we note that outside of the 
relativistic regimes, the HCL and 75 lattice behave very 
similarly as they both feature alternation of non-trivial 
and trivial phases as a function of E^. This important 
result emphasizes that the optical lattice topology dra- 
matically influences the nature of the insulating phases 
within the relativistic regime. Consequently, engineered 
gauge fields offer the unique possibility to explore Z 2 - 
phase transitions as the optical-lattice topology is modi- 
fied. 

Different schemes are available for distinguishing the 
QSH and normal insulating phases in a cold- atom ex- 
periment [12]. In absence of spin-mixing perturbations 
the spin is conserved and the Z2 topological constant 
reduces to v — |A^|mod2 because the individual spin- 
components are associated to N^ = —N±. Therefore v 
can be determined through density measurements and 
by applying the Streda formula [31]. Alternatively, the 
presence of helical edge-states can be also probed through 



light-scattering methods by exploiting the method re- 
ported in Ref. [32]. 



V. CONCLUSIONS 

We have studied the effect of Haldane-type SOI on 
two interrelated two-dimensional lattices: the honey- 
comb and the Tz lattices. Both can be transformed into 
each other by selectively switching on/off specific lattice 
sites. In absence of SOI, the main effect of the transition 
from HCL to 73 is the appearance of a localized band 
at the charge neutrality point. For finite SOI though, a 
topological phase transition occurs that corresponds to 
the destruction of the QSH phase which is present in the 
HCL but not allowed in the T3 . Since the occurrence 
of the QSH phase is related to the lattice topology, it 
is already destroyed by adding an infinitesimal weakly 
coupled lattice site to the HCL which gives rise to the 
7i lattice. Finally, we have proposed a method for im- 
plementing this topological phase transition for fermionic 
cold atoms trapped in optical lattices. In this context, 
the QSH phases are obtained by synthesizing a gauge 
field reproducing the effects of the Haldane-type SOI for 
electrons. 
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Appendix A: The bulk energy spectrum 

The energy spectrum of the 73 lattice with Haldane- 
type [18] SOI is obtained by diagonalization of the tight- 
binding Hamiltonian % = Ho +^so- Here the first term 
is a nearest-neighbor hopping tight-binding Hamiltonian 
and includes the spin independent nearest-neighbor hop- 
ping amplitude t. The second term — defined in Eq. (1) 
- describes the Haldane-type SOI. Because of transla- 
tional invariance, this Hamiltonian can be diagonalized 
in reciprocal space. The evaluation of the contribution 
due to the Haldane-SOI requires some care because of 
the different connectivity of hubs and rims. A second- 
next neighbor hopping process always interconnects sites 
of the same species (either rim A, rim B or hub). Be- 
cause of the two possible distinct processes that connect 
two given hubs, the contributions arising from these two 
paths cancel. Therefore, the Haldane-SOI effectively only 
couples the rims of each species but not the hubs. 



The Hamiltonian now can be written in matrix form 

/0 A* N 
H[k]=t\A A 
\0 A* 



where we have defined the coefficients 




_ifc-vi 



Jk-v 2 



A=l + e" v v 1 + e" 

B = 2{sin(fc • vi) — sin(fe • V2) 



■sin[fe- (vi +v 2 )]}. 



Here the translational vectors of the T3 lattice are de- 
fined as vi = (3/2; -\/3/2)4 and v 2 = (3/2; \/3/2)4. 
Equation (Al) corresponds to a six component spinorial 
wave function ^fe(r) = (^a, ^h? ^b) ® cr, where we have 
introduced the components of the wave function on the 
two rims and the hub sites and the spin contribution a. 
Diagonalization of H[k] yields the energy spectrum for 
the 75 lattice in presence of Haldane-SOI 

e (k) =0 (A2a) 

e±(k) = ± [2 [3(t 2 + tl ) + 2t 2 (cos[k • vi] (A2b) 

+ cos[k • (vi - v 2 )] + cos[k • v 2 ]) 
— ts Q (cos[2k • vi] — 2 cos[k • v 2 ] 
+ cos[2fc • v 2 ] + cos[2fc • (vi — v 2 )] 
+ 2 cos[2k • vi - k • v 2 ] + 4(sin[fc • vi] 

>, 1/2 
+ sin[fe • (vi - v 2 )] sin[fc • v 2 ]))] > 

where k is the two-dimensional momentum vector and 
each band is two- fold degenerate. 

Appendix B: Long-wavelength approximation 

For momenta close to the two independent if-points 

— K = 2tt/3£ (1,-- n/3/3) and K' = 2tt/34(1,+>/3/3) 

— we can apply a long-wavelength approximation to 
the Schrodinger equation underlying the Hamiltonian 
^0 + ^so- This approximation consists in expressing 
the spatial part of the wave function as the product of a 
fast- varying part times a slow- varying part. In absence of 
perturbations that induce a mixing of the two K points, 
we consider a single if -point — namely K — and the 
wave function can be written as 



* a ( J R a )oce ifeR °Ff( J R a ). 



(Bl) 



Here a G {A,B,H} and R a is the lattice site coordinate. 
We substitute this wave-function into the Schrodinger 
equation and use the expansion 

FfiRc^dj) * Ff(R a ,)±d r V r Ff(r)\ r=R +0(\d\ 2 ) ■ 

This approximation is valid for energies close to the 
charge neutrality point where the system energy is al- 
most zero. In this case the associated Fermi wavelength 



Ap ~ l/|e| is bigger than the modulus of the displace- 
ment vectors Ap ^> \dj\. Collecting all the terms we 
are left with the expression (2) where vp = 3£ot/2 and 
Aso = V3tso. 

The pseudo-spin matrices T, have been defined in (3). 
These matrices fulfill the algebra of the angular momen- 
tum [E^,Ej] = ieijk^k and form a 3-dimensional repre- 
sentation of SU(2). However, contrary to the Pauli matri- 
ces, they do not form a Clifford algebra, i.e., {£$, Sj} ^ 
2SijIs. By introducing a rotation operator around the z 
axis defined by T> z {(j)) = exp (— iE^0), a generic state \a) 
is transformed into itself by T> z (27r)\a) —> \a), implying 
that the pseudo-spin X) describes an integer spin 5=1. 




(b) 




FIG. 5: (Color online). Peierls phases generating the SU(2) 
gauge field for (a) the honeycomb lattice and (b) in the 75 lat- 
tice. 



Appendix C: Topological properties of the 75 lattice 
with Haldane-SOI 



The Z2 topological invariant characterizes the change 
in the Kramers degeneracy of a time-reversal symmet- 
ric system [2]. It can be deduced from the topologi- 
cal structure of the Bloch wave functions of the bulk 
crystal in the first Brillouin zone. Here, there are spe- 
cial points k = Mi that are time-reversal invariant and 
satisfy — M^ = Mj + G for a reciprocal-lattice vector 
G. For the Ts lattice we have four of these points. 
These are defined by M i=(ni5n2 ) = l/2(n 1 g 1 + n 2 g2) 
with rii, 77,2 G {0,±1}, where we have introduced the 
reciprocal lattice vectors gi = 2tt/£o( 1/3; — a/3/3) and 
g2 = 2tt/£o (1/3; + >/3/3). Defining the time-reversal 
symmetry operator as T = ex.p(i7ra y /2)JC with the op- 
erator of complex conjugation /C, the four M points are 
found to fulfill H(Mi) = rH(M i )r~ 1 . 

In the case of the 75 lattice the system is invariant un- 
der parity operation, which corresponds to an exchange 
of rim A with rim B and is expressed by the operator 



V 




(CI) 



This implies that [H(Mi),V] = for all four time-reversal 
invariant momenta. According to Fu and Kane [23] the 
Z2 topological invariant v is related to the parity eigen- 
values ^2m(M^) of the 27Ti-th occupied energy bands at 
the four time-reversal invariant points M^ through the 
relation 



N 



(-ir=nn wm*). 



(C2) 



We have evaluated the eigenstates of H in the points M^ 
analytically and therefore determined straightforwardly 
the parity eigenvalues of the occupied bands. At filling 
1/3 and 2/3 we find that all four eigenvalues are positive 
determining therefore a trivial Z2 number v = 0. 



Appendix D: Subjecting the Honeycomb and 
75 lattices to an SU(2) gauge field 

On a lattice, the gauge field A enters the formalism 
through the Peierls (or Berry's) phases: as a particle 
hops from a site i to a nearest-neighboring site j, the 
wave function acquires a non-trivial phase e l6> ^, where 
6ij — J? A • dZ. The set of Peierls phases {%}, defined 
on the whole lattice, therefore witnesses the presence of 
a specific gauge field A. The physical observables asso- 
ciated to the gauge field, or equivalently to the set {#r/}, 
are given by the loop operators defined on each plaquette 
□ : 



W(B) 



J9 12 „i9 2 3 



(Dl) 



where the plaquette □ is delimited by the sites 
{1,2,..., N}. In the case of an Abelian or U(l) gauge 
field, the loop operator is related to the magnetic flux 
$(□) penetrating the plaquette through the relation 

In this work, we considered a gauge field of the follow- 
ing form 



A = - (0, 27T(j)X, 0) 

o 



)(7 Z = 



A^ 








(D2) 



that does not couple the two spin species together. Hence 
up and down spins feel individual and opposite magnetic 
fluxes 



* t (D) = / V 

= 2tt(/>, 



■ / V X A; 

Ju 



-^(D) (D3) 



that are uniform. The loop operators corresponding to 
the gauge field (D2) thus yield W = e l27V(f)az within every 
plaquette of the lattice. The helical edge-states, hallmark 
signature of the QSH phase, are a direct consequence of 
the double-Hall system produced by the specific gauge 
field (D2) which preserves time-reversal symmetry. 



The gauge field (D2) can be engineered in an optical- 
lattice experiment by inducing the appropriate set of 
Peierls phases {0^-} such that 



W 



Q W 12^1023 



, . e 



(D4) 



for all the plaquettes constituting the honeycomb and 
73 lattices. Note that the Peierls phases can be generated 
through different schemes, based on the Raman transi- 
tions produced by additional lasers or fields along the 
appropriate links [12, 33, 34]. An adequate set of Peierls 
phases {0ij} satisfying Eq. (D4), and thus reproducing 
the effects of the gauge field (D2) on the honeycomb and 
T3 lattices, are represented in Fig. 5. 

In the case of the HCL, a single space-dependent 
Peierls phase Q(m) = 2i\(\)Vfia z is necessary, as illustrated 
in Fig. 5a. The single-particle Schrodinger equation then 
takes the form of two coupled Harper difference equa- 
tions for the wave-function \I/ = (^a^b) defined at the 
site (m, n), 

^a(^i, n) =?/ ; b(^ n — 1) + i)j&{m + l,n — 1) (D5a) 

+ e i2 ^ ma ^ B (m,ri), 
V^b(^, n) =^a(^, n H- 1) + ^a(^ — 1, n + 1) (D5b) 

+ e- i2 ^ m ^A(m,n), 



and the loop operator indeed equals 
within all the hexagonal plaquettes. 



W 



.i27vd>a z 



In the case of the T3 lattice, four space-dependent 
Peierls phases are required, 

01 (ra) = 27T0(3ra - l/2)a z , (D6a) 

<9 2 (m) = 27r0(3m - l)a z , (D6b) 

6 3 (m) = 27r0(-3m - l/2)a z , (D6c) 

4 (m) = 27T0(-3ra - l)a z , (D6d) 



as represented in Fig. 5b. The Harper equations for the 
wave-function ^ = (^a, ^b, ^h) defined at the site (m, n) 
yields 



i/>A(m, n) =e"^ 3(m VH(ra, n ) + ^ H ( m , n - 1) (D7a) 

+ e-^ 2(m+1 VH(m + l,n-l), 
ip B (m, n) =e _i0l(m VH(m, n) + ?MO, n + 1) (D7b) 

+ e- i ^ m - 1 ^(m- l,n + l), 
^ H (m, n) =e i03(m VA(m, n) + ^(m, n - 1) (D7c) 

+ e i0l(m VB(^i, ™) + ?AaO, n + 1) 

+ e^ 4(m VB(m + l,n-l) 

+ e^ 2(m VA(m-l,n + l) 



and one verifies that W(U) = e i27T(f)az 
bus forming the T3 lattice. 



within each rhom- 



The energy spectra of the honeycomb and T3 lat- 
tices are readily obtained by diagonalizing the Harper 
equations (D5) and (D7) numerically, leading to the 
butterfly-like structures represented in Figs. 4a and b 
of the manuscript. Since the gauge field does not cou- 
ple the two spin components, the Z2-topological invari- 
ant v — I iV^ I mod 2 is directly obtained by computing 
the Chern number iV-f associated to the up-spin. This 
Chern number can be evaluated numerically by consid- 
ering one copy of the two-fold Harper equations [30]. 
Considering the more general relation v = n a mod2, 
the Z2 -topological invariant can equally be evaluated by 
computing the spin-Chern number n a associated to the 
two-fold Harper equations through the twisted-periodic- 
boundary method [21]. 



[1] K. v. Klitzing, G. Dorda, and M. Pepper, Phys. Rev. 

Lett. 45, 494 (1980). 
[2] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 146802 

(2005); Phys. Rev. Lett. 95, 226801 (2005). 
[3] B. A. Bernevig and S.-C. Zhang, Phys. Rev. Lett. 96, 

106802 (2006). 
[4] B. A. Bernevig, T. L. Hughes and S.-C. Zhang, Science 

314, 1757 (2006). 
[5] M. Koenig et a/., Science 318, 766 (2007). 
[6] A. Roth et al, Science 325, 294 (2009). 
[7] Y. Xia et a/., Nature Phys. 5, 398 (2009). 
[8] H. Zhang et aL, Nature Phys. 5, 438 (2009). 
[9] C. Xu and J. Moore, Phys. Rev. B. 73, 045322 (2006). 
[10] I. Bloch, J. Dalibard and W. Zwerger, Rev. Mod. Phys. 

80, 885 (2008). 
[11] M. Lewenstein et a/., Adv. Phys. 56, 243 (2007). 
[12] N. Goldman et a/., Phys. Rev. Lett. 105, 255302 (2010). 
[13] T. D. Stanescu et a/., Phys. Rev. A 82, 013608 (2010). 



[14] L. Sanchez-Palencia and M. Lewenstein, Nature Phys. 6, 

87 (2010). 
[15] B. Sutherland, Phys. Rev. B 34, 5208 (1986). 
[16] J. Vidal, R. Mosseri, and B. Doucot, Phys. Rev. Lett. 

81, 5888 (1998); J. Vidal, P. Butaud, B. Doucot, and 

R. Mosseri, Phys. Rev. B 64, 155306 (2001). 
[17] S. E. Korshunov, Phys. Rev. B 63, 134503 (2001); 

V. Cataudella and R. Fazio, EPL 61,341 (2003). 
[18] F. D. H. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 
[19] D. Bercioux, D. F. Urban, H. Grabert and W. Hausler, 

Phys. Rev. A 80, 063603 (2009). 
[20] D. Bercioux, M. Governale, V. Cataudella and V. M. Ra- 

maglia, Phys. Rev. Lett. 93, 056802 (2004); Phys. Rev. 

B 72, 075305 (2005). 
[21] D. N. Sheng, Z. Y. Weng, L. Sheng and F. D. Haldane, 

Phys. Rev. Lett. 97, 036808 (2006). 
[22] T. Fukui and Y. Hatsugai, Phys. Rev. B 75, 121403(R) 

(2007). 



[23] L. Fu and C. L. Kane, Phys. Rev. B 76, 045302 (2007). 205414 (2006). 

[24] C. Weeks and M. Franz, Phys. Rev. B 82, 085310 (2010). [31] R. O. Umucalilar, H. Zhai and M. 6. Oktel, Phys. Rev. 

[25] H.-M. Guo and M. Franz, Phys. Rev. B 80, 113102 Lett. 100, 070402 (2008). 

(2009). [32] X.-J. Liu, C. Wu and J. Sinova, Phys. Rev. A 81, 033622 
[26] W. Li and R. Tao, arXiv:1001.4168, unpublished. (2010). 

[27] K. Lee et a/, Phys. Rev. A 80, 043411 (2009). [33] D. Jaksch and P. Zoller, New J. Phys. 5, 56 (2003). 

[28] M. Rizzi, V. Cataudella and R. Fazio, Phys. Rev. B 73, [34] K. Osterloh, M. Baig, L. Santos, P. Zoller, and M. Lewen- 

144511 (2006). stein, Phys. Rev. Lett. 95, 010403 (2005). 

[29] A. Bermudez et a/., New J. Phys. 12, 033041 (2010). 
[30] Y. Hatsugai, T. Fukui and H. Aoki, Phys. Rev. B 74, 



